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We investigate the formation of localized structures with a varying width in one and two- 
dimensional systems. The mechanism of stabilization is attributed to strong nonlocal coupling 
mediated by a Lorentzian type of Kernel. We show that, in addition to stable dips found re¬ 
cently [see, e.g., C. Fernandez-Oto, M. G. Clerc, D. Escaff, and M. Tlidi, Phys. Rev. Lett. 110 , 
174101 (2013)1, exist stable localized peaks which appear as a result of strong nonlocal coupling, 
i.e. mediated by a coupling that decays with the distance slower than an exponential. We ap¬ 
plied this mechanism to arid ecosystems by considering a prototype model of a Nagumo type. In 
one-dimension, we study the front that connects the stable uniformly vegetated state with the bare 
one under the effect of strong nonlocal coupling. We show that strong nonlocal coupling stabilizes 
both—dip and peak—localized structures. We show analytically and numerically that the width 
of localized dip, which we interpret as fairy circle. Increases strongly with the aridity parameter. 
This prediction is in agreement with filed observations. In addition, we predict that the width of 
localized patch decreases with the degree of aridity. Numerical results are in close agreement with 
analytical predictions. 

PACS numbers: 


I. INTRODUCTION 

Localized structures (LS’s) in dissipative media have been observed in various field of nonlinear science such as 
fluid dynamics, optics, laser physics, chemistry, and plant ecology (see recent overviews [1, 2]). Localized structures 
consist of isolated or randomly distributed spots surrounded by regions in the uniform state. They may consist of 
dips embedded in the homogeneous background. They are often called spatial solitons, dissipative solitons, localized 
patterns, cavity solitons, or auto-solitons depending on the physical contexts in which their were observed. Localized 
structures can occur either in the presence [3] or in the absence [4] of a symmetry breaking instability. In the last 
case, bistability between uniform states is a prerequisite condition for LS’s formation. However, in the presence of 
symmetry breaking instability, the coexistence between a single uniform solution and a patterned state allows for the 
stabilization of LS’s [4]. In this case, bistability condition between uniform solutions is not a necessary condition for 
generating LS’s [5]. 

Spatial coupling in many spatially extended systems is nonlocal. The Kernel function that characterizes the non¬ 
locality can be either weak and strong. If the Kernel function decays asymptotically to infinity slower than an 
exponential function, the nonlocal coupling is said to be strong [6]. If the Kernel function decays asymptotically to 
infinity faster than an exponential function, the nonlocal coupling is said to be weak. Self-organization phenomenon 
leading to the formation of either extended or localized patterns under a local and nonlocal coupling occur in various 
systems such as fluid dynamics [7], firing of cells [8, 9], propagation of infectious diseases [10], chemical reactions 
[11, 12], population dynamics [13-15], nonlinear optics [16-22], granular [23], and vegetation patterns [24-28]. 

We focus on bistable regime far from any symmetry breaking instability, i.e. far from any Turing instability. In 
this case, the behavior of many systems is governed by front dynamics or domains between the homogeneous steady 
states. When the nonlocal coupling is weak, the interaction between two fronts is usually described by the behavior 
of the tail of fronts. However, for strong nonlocal coupling, the interaction is controlled by the whole Kernel function 
and not by the asymptotic behavior of the front tails [6]. 

When a non-local coupling is weak, the asymptotic behavior of front solutions is characterized by exponential decay 
or damping oscillations. In the former case, front interaction is always attractive and decays exponentially with the 
distance between the fronts. In two dimensional setting, LS’s resulting from fronts interaction are unstable. In the 
case of damping oscillations, front interactions alternate between attractive and repulsive with an intensity that decays 
exponentially with the distance between the two fronts [29]. For a fixed value of parameters, a family of stable one 
dimensional localized structures with different sizes has been reported in [13, 30-32]. An important difference appears 
when considering a strong nonlocal coupling. This difference is that the interaction between fronts can be repulsive 
[33, 34]. 

Strong nonlocal coupling has been observed experimentally in various systems. Indeed, several experimental mea- 
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Figure 1: (Color online) Localized vegetation patterns: (a) An aerial photo taken on 05 March 2010 shows fairy circles in 
the Namibrand region, Namibia. Image courtesy of Professor Norbert Juergens. (b) An aerial photograph showing localized 
vegetation patches, Zambia [61]. 


surements of nonlocal response in the form of Lorentzian or a generalized Lorentzian have been carried out in nematic 
liquid crystals cells [36, 37]. Experimental reconstruction of strong nonlocal coupling has also been performed in 
photorefractive materials [38]. In this case, the strong nonlocal coupling is originating from the thermal medium 
effects. In population dynamics such as vegetation, it has been shown experimentally that seed dispersion may be 
described by a Lorentzian [39] . 

We consider a proto model for population dynamics, namely the strong nonlocal Nagumo equation [33]. This 
model possesses two relevant properties: strong nonlocal coupling and bistability between uniformly vegetated and 
bare states. We focus on a regime far from any symmetry breaking or Turing type of instabilities. In this regime 
localized structures resulting from strong nonlocal coupling can be stabilized in a wide range of parameters [33]. 
We will investigate two types of localized vegetation structures: i) isolated or randomly distributed circular areas 
devoid of any vegetation, often called fairy circles (EC’s), and ii) isolated or randomly distributed circular areas of 
vegetation, surrounded by a bare region. An example of EC’s is shown in an areal photograph (see Fig. la). They 
are observed in vast territories in southern Angola, Namibia, and South Africa [40-48]. The size of these circles 
can reach diameters of up to 12m. An in-depth investigation of several hypotheses concerning their origin has been 
performed by van Rooyen et al. [49]. In this study, these authors have been able to excluded the possible existence 
radioactive areas inapt for the development of plants, the termites activity and the release of allelopathic compounds. 
We attribute two main ingredients to their stabilization: the bistability between the bare state and the uniformly 
vegetation state, and Lorentzian-like non-local coupling that models the competition between plants [34] . We provide 
detailed analysis of the fairy circle formation in a simple population dynamics model. In addition we show that 
the above mechanism applies to another type of localized structures that consist of isolated or randomly distributed 
vegetation patches surrounded by a bare state. An example of this behavior is shown in Fig. lb. In this paper, we 
investigate analytically and numerically the formation of both—fairy circles and localized patches-and their existence 
range. Our theoretical analysis shows that exist a Maxwell point above which localized patches are stable, while below 
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this point fairy circles appear. Finally, we investigate how the degree of aridity affects the width of both types of 
localized vegetation structures. 

This paper is organized as follows. After a briefly introducing the model describing the vegetation dynamics, 
namely the Nagumo model with strong nonlocal coupling mediated by a Lorentzian function (Sec. II). We describe 
the dynamics of a single front in one dimension, and its asymptotic behaviors (Sec. III). The analytical and numerical 
analysis of the interaction between fronts connecting the uniformly vegetated and the bare steady states is described 
in Sec. IV, where we discuss the formation of both fairy circles and localized patches. Close to the Maxwell point, we 
derive a formula for the width of both localized structures as a function of the degree of the aridity in one dimension. 
We conclude in Sec. V. 


II. THE NAGUMO MODEL 


Several models describing vegetation patterns and self-organization in arid and semiarid landscapes have been 
proposed during last two decades. They can be classified into three types. The first approach is based on the 
relationship between the structure of individual plants and the facilitation-competition interactions existing within 
plant communities [24-26, 50, 51]. The second is based on the reaction-diffusion approach which takes into account of 
the influence of water transport by below ground diffusion and/or above ground run-off [52-58]. The third approach 
focuses on the role of environmental randomness as a source of noise induced symmetry breaking transitions [59-62] . 
Recently, the reduction of a generic interaction-redistribution model, which belong to the first class of ecological 
type of models [26], to a Nagumo-type model has been established [34] . Here, we consider the variational nonlocal 
Nagumo-type equation 


dtu 


u{a — u){u — 1)+ V^u + eu / u^(r-|-r', t)Ar(r')(ir' 

Jn 


( 1 ) 


where u (r, t) is a normalized scalar held that represents the population density or biomass, a is a parameter describing 
the environment adversity or the degree of aridity, t is time. We consider that population or plant community 
established on a spatially uniform territory fl. The vegetation spatial propagation, via seed dispersion and/or other 
natural mechanisms, is usually modeled by a non-local coupling with a Gaussian like Kernel [24-26]. For simplicity, we 
consider only the first term in Taylor expansion of the dispersion. This approximation leads to the Laplace operator, 
= dxx + dyy acting in the space r = {x,y). The last term in Eq. 1 describes the competitive interaction between 
individual plants through their roots. The nonlocal coupling intensity, denoted by e, should be positive to ensure a 
competitive interaction between plants. The kernel function has the form K (r) = 6 (r) — (r), where 6 (r) is the 

delta function. The inclusion of the S function in the analysis allows us to avoid the variation of the spatially uniform 
states as a function the nonlocal intensity e, and 


fa (r) 


Nn 

i + (M/a)"’ 


( 2 ) 


which has an effective range a. For the sake of simplicity, we consider that the length cr is a constant, independent 
of the biomass. Hence, we assume all plants have the same root size, that is, we neglect allometric effect [50, 63] . At 
large distance, the asymptotic behavior of the Kernel is determined by n, and Nn is normalization constant. 

It is worth to emphasize that, since the strong nonlocal term in Eq. (1) models the interaction between individual 
plants at the community level , we refer to it as strong nonlocal interaction or coupling. In contrast, non-localities 
describing transport processes, as seed dispersion [24-26], for sake of simplicity, we are modeling by the Laplace 
operator. 

The Eq. (1) is variational and it is described by 


dtU 


5F dF „ 


( 3 ) 


where is a Lyaponov functional that can only decrease in the course of time. Accordingly, any initial distribution 
u (r, t) evolves towards a homogeneous or inhomogeneous (periodic or localized) state corresponding to a local or 
global minimum of F. The Lyaponov functional reads 


e 



■^M= y'^|^(|Vu|Vu(u)|dr 

f f vf{v)u^{v')K{r — v')drdT' 

J qJ Q 
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and 

V{u) = ^{u- 1)^ + y (a - 1/2) (3 - 2u) . (4) 

The Eq. (1) admits three spatially uniform solutions u = 0, u = a and u = 1. The bare state u = 0 is always 
stable, and represents not plant state. The uniform state u = a is always unstable. For large values of a, the climate 
becomes more and more arid. The uniformly vegetated state u = 1 may undergo a may undergo a symmetry breaking 
type of instability (often called Turing instability), that leads to pattern formation. In one dimensional system, and 
for n = 2, the threshold for that instability satisfies 

P = ea^exp{-l3), ( 5 ) 

with P = ^J\ + cr^(2e + a — 1) — 1. In what follows, we will focus on a regime far from any pattern forming instability. 


III. FRONTS 

We consider a bistable regime where u = 0 and u = 1 are both linearly stable. From Eq. (4), F (0) =0 and 
V (1) = (a — 1/2) /6. Therefore, when a < 1/2, the most favorable state is the uniformly vegetated one. When 
a > 1/2, the bare state is more stable than the uniformly vegetated one. There exist a particular point where both 
states are equally stable. This point is usually called the Maxwell point [35], and corresponds to a = 1/2. 

Depending on the value of a, front connecting both state will propagate towards the most stable state. An 
example of a single front is illustrated in Fig. 2. a. The time-space diagram of Fig. 2.b shows how the most stable 
state, corresponding to the bare state (a > 1/2), invades the uniformly vegetated state, with a constant speed. Fronts 
propagate following the minimization of the potential (4) and the front velocity is proportional to the energy difference 
between equilibria V (1) — V (0). The front is motionless at the Maxwell point for a = 1/2. At this point and in the 
absence of nonlocal coupling, e = 0, front solutions read 

u± (x -xo) = ^ tanh > (6) 

where xq corresponds to the interphase position. The front links the barren state from x = —oo to the uniformly 
vegetated state at a: = oo. The opposite connection correspond to u_. 

The asymptotic behavior of the front solutions (6) obeys to an exponential law in the form 

u+ (x « Xo) « 

u+ (x » a;o) « 1 - 

Around the barren state (u = 0), the inclusion of the nonlocal interaction does not modify the asymptotic behavior 
of the front, since the nonlocality in model (1) is nonlinear. Let us examine the effect of nonlocal term around the 
uniformly vegetated steady state (u = 1). For this purpose, let us assume that the asymptotic behavior of the front 
obeys to an exponential law in the form 

u (a: » Xo) « 1 - 

where c is a constant, and the exponent A obeys the equation 

— - -I- e (3 — 2gcr (A)) = 0, 


where 

/ OO 

cosh (Ax) (x) dx. 

-OO 

This equation has been successfully used to explain the emergence of localized domains. When A-solutions have 
non-null imaginary part, spatially damped oscillations on the front profile are induced, leading to the stabilization of 
localized domain. This mechanism is well documented for either local or nonlocal systems [13, 30, 31, 64[. 

The function exists when the Kernel decays faster than an exponential one, i.e., a weak nonlocal coupling. In 
the case of strong nonlocal interaction, diverges and the above analysis is no more longer valid. 
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Figure 2: (Color online) Front propagation obtained from numerical simulations of Eq. (1). (a) Biomass front profile, (b) 
Spatiotemporal diagram that shows the front movement with a constant speed. The parameters are a = 0.51, a = 0.7, n = 2 
and e = 1.0. 


To determine the asymptotic behavior of the front around the uniformly vegetated state under strong nonlocal 
coupling, we perform a regular perturbation analysis in terms of small parameter e. At the Maxwell point, we expand 
the field as 


u {x) = uq {x) + eui {x) + e^U 2 {x) + ..., (7) 

where uq = is the motionless front provided by the Eq. (6). Replacing (7) in Eq.(l), and making an expansion in 
series of e. At order e, we obtain: 


{dxx - 1/2 + 3uo (1 - uo)} ui = 

{uo ul{x + x')fa{x')dx' - ug) . 

Let us focus in the region x ^ Xq, and neglecting all exponential corrections coming from uq, then we obtain 

nOO 

{dxx - 1/2} ui= - {x') dx'. (8) 

J X — Xq 

Equation (8) is a linear inhomogeneous equation for the correction ui. For strong nonlocal coupling, the particular 
solution of Eq. 8 dominates over the homogeneous one which is exponentially small. For instance, if we consider a 
Kernel like 

N 

fa {x) « — for x » 1, (9) 

x” 

with n > 1 and N a normalization constant. For x — Xq ^ 1, the front approaches asymptotically the following 
solution 


2eN 

{n - 1) (x - Xo)”“^ 


( 10 ) 


This solution decays according to a power law 1 —n. To check the power obtained from the above analysis, we perform 
a numerical simulations of the full model Eq. (1). The result obtained from Eq. (9) and the numerical simulations 
are shown in in Fig. 3. Both results agree perfectly without any adjusting parameter. Note that both analytical 
calculations and numerical simulations predict the existence of one peak in the spatial profile of the front. This peak 
takes place at the interface separating both homogeneous steady states as shown in Fig. 3. 
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Figure 3: (Color online) Numerical solutions of model (1), with the Lorentzian kernel (2). (a) Front profile at the Maxwell 
point a = 0.5, other parameters are a = 2, n = 2 and e = 0.35. (b) Decay of front to the uniformly vegetated state for a = 0.5, 
a = 2, n = 2 and e = 0.1. In (b), the upper panel shows the comparison between the numerical calculation (dotted line) and 
the analytical estimation (10) (continuous line), while the lower panel displays the same comparison in a log-log plot. 


IV. LOCALIZED VEGETATION PATTERNS 

Far from a symmetry breaking instability, localized structures can be stable as a results of front interactions. This 
phenomenon occurs when the spatial profile of the front exhibits damped oscillations [13, 64]. However, around 
the bare state damped oscillations are nonphysical since the biomass is a positive defined quantity. A stabilization 
mechanism of LS’s based on combined influence of strong nonlocal coupling and bistability has been proposed [33]. 
This mechanism has been applied to explain the origin of the fairy circles phenomenon in realistic ecological model 
[34]. To the best our knowledge, there is no other analytical understanding of the circular shape of the fairy circle. 
We have shown in addition that the diameter of the single fairy circle is intrinsic to the dynamics of the system such 
as the competitive interaction between plant and the redistribution of resources [34] . We believe that extrinsic causes 
such as termite or ant or others external environment cannot explain the circular shape of fairy circles. 

In this section we provide a detailed analysis of front interactions leading to stabilization of both fairy circles and 
localized patches. For both types of localized vegetation structures, we discuss how the level of aridity affects their 
diameter. 
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Figure 4: (Color online) Stationary fairy circle obtained from numerical simulations of Eq. (1). (a) One dimension with 
a = 0.492, a — 0.8, n = 2 and e = 0.5. (b) Two dimensions with a = 0.38, a = 0.2, n = 2.2 and e = 1. 


A. Fairy circles 

The model Eq. (1) admits stable localized structures in the form of bare state embedded in herbaceous vegetation 
matrix. An example of such behavior is shows in Fig. 4. They are stable and permanent structures. A single fairy 
circle exhibits a fringe formed by tall grass that separate the bare state to the uniformly vegetated as shown in Fig. la. 
From numerical simulations, we see that the biomass possess one peak that takes place in between the bare and the 
uniformly vegetated states as shows in Fig. 4. This can be explain by the fact that inside the circle the competition 
between plant is small. Indeed the length of the plant root is much smaller than the diameter of the fairy circle. 

We analytically investigate the formation of a single fairy circle in one spatial dimension. We focus on the parameter 
region near to Maxwell point (^ — a = ry <C 1) and we consider a small non-local coupling (e <C 1). We look for a 
solution of Eq. (1) that has the form of a slightly perturbed linear superposition of two fronts 

u{x, t) = u+{x — S{t)) -I- u_ (a; -I- 6(t)) + W, (11) 

where W(x,u^,U-), dtS, tj and e are small in a sense that will be pinned down bellow. While u± are defined by Eq. 

( 6 ). 

Replacing anstatz (11) in Eq. (1) and neglecting high order terms in e, we obtain 

— {dxU^ — dxU-)dtS = u+u_(3 — 3it+ — 3 m_) 

/ OO 

v?‘{x + x', t)K{x')dx 

-OO 

—r]{u+ + U- — l)(u+ + U-) -b LW -b h.o.t., (12) 

where the linear operator L has the form 

L = — ^ -b 3u+ — 3u^ -b dxx + 3m- — 3u?_ — 6u+u_. 

To solve the above equation, we consider the inner product {g\h) = {g{x)h(x))dx. Then, the linear operator L_|_ 

is self-adjoint and LdxU± ~ 0. The solvability condition gives 


{dxu+\dxu+) dt5 = G + rj {dxu+\{l - u+)u+) 


( 13 ) 
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Figure 5: (Color online) Width of a fairy circle, as a function of aridity a in one dimension. The solid line represents the 
analytical prediction from Eq. (16). The squares are obtained by numerical simulations of Eq. (1). The inset is the spatial 
proHle of biomass. Parameters are ct = 0.8, n = 2 and e = 0.5. 


where 


G{u±,a) 


dxU+\u+ 



u^{x', t)K{x' 



We neglect also the terms smaller than 1/(5^" ^ in Eq. (13). Then, it is obtained 


‘ (n-l)(2<5)"-i 


V2r), 


(14) 


Strictly speaking, this equation of motion is quantitatively valid when dt6 ^ ^ 77 <?; 1. 

This result is valid for any power n > 1 in the Kernel function (2). We consider n = 2 and iV„ = I/ttct in one 
dimension. Then, the equation (14) reads 



where A = 25 is the width of the LS. The equilibrium width is given by 

* 3ecr _ 

Aeg = jj Y' (16) 

^(2 -«) 

The linear stability analysis allows us to determine the eigenvalue A = —2^717]'^ jSea. Therefore, for competitive 
interaction, i.e., e > 0, fairy circles are always stable. 

We plot the width of the FC in one dimension as function of the aridity a in Figure 5. The width of localized 
structures grows as aridity increases. At the Maxwell point, i.e., a = 1/2, the width of the fairy circle becomes infinite. 
In order to check the approximations used to derive the equilibrium size (16), we perform numerical simulations of 
Eq. (1). Both results are in good agreement, without any fitting parameter. 

The fairy circles are observed in vast territories in southern Angola, Namibia, and South Africa [41, 42], where the 
annual rainfall ranges between 50 and 150 mm [49]. The size of FC increases from South to North where the climate 
becomes more and more arid [49] . The size can also be affected by the rainfall and nutrients [43]. Fairy circles average 
diameters varies in the range of 2 m-12 m [49]. In agreement with field observations. Fig. 5 shows indeed that the 
fairy circles diameter increases with the aridity. 

Therefore, one-dimensional front interaction explains why the fairy circles size increase with the aridity. To wit, as 
environment aridity increases, the bare state becomes more and more favorable, increasing the fairy circles size. This 
mechanism demands, however, that the uniformly vegetated state must be always the most favorable one (a < 1/2), 
otherwise the bare state propagates indefinitely. The same tendency is observed in two-dimensional simulations [34] . 
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B. Localized vegetation patch 

In this subsection, we investigate the formation of a single localized patch that consist of a circular vegetated state 
surrounded by a bare state. This behavior occurs for high values of the aridity parameter, i.e., a > 1/2. An example 
of a single localized patch is illustrated in Fig. 6. This localized solution corresponds to the counterpart of fairy 
circles. 

Following a similar strategy than in the previous subsection (small e and rjf = a — 1/2), the patch width A in one 
spatial dimension obeys 


6V2eA„cr” 
(n- l)(A)"’-i 


2\/2r]f. 


(17) 


Note that, like in Eq. (14), the result obtain in Eq. (17) is generic for any n in (2). For n = 2, the stable vegetated 
patch have the size 

A 3ecr 

Aeg = p-. (18) 

2) 

The formula (18) is plotted in Fig. 7 by a solid line. Confrontation with direct numerical computation for the 
localized patch width is in good agreement, as shown in Fig. 7. There is no available data from the field observation 
to confirm this theoretical prediction. 


C. Bifurcation diagram 

In this subsection, we establish the bifurcation diagram for both types of localized vegetation structures. We fix 
the length of the completion between plant a, and we vary the degree of aridity a and the strength of the competitive 
interaction e. We numerically establish a stability range of a single fairy circle and the localized patch in one 
dimensional setting. This analysis is summarized in the parameter plane {a,e) of Fig. 8. For a < 1/2, a single fairy 
circle is stable in the region FC as indicated in Fig. 8. This stability region is bounded from below by dots and bounded 
from left by the Maxwell point (a = 1/2). Dynamically speaking, dots correspond to a saddle-node bifurcation. The 
parameter zone A indicates the regime where a FC shrinks and disappears. For large values of the strength of the 



Figure 6: (Color online) Stationary localized patch obtained from numerical simulations of Eq. (1). (a) One dimension with 
a = 0.52, a = 0.8, n = 2 and e = 0.5. (b) Two dimensions with a = 0.56, a = 0.1, n = 2.5 and e = 10. 
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Figure 7: (Color online) Width of a localized vegetation patch, as a function of aridity a in one dimension. The solid line 
represents the analytical prediction from Eq. (18). The dots are obtained by numerical simulations of Eq. (1). The inset is the 
spatial profile of biomass. Parameters are o = 0.8, n = 2 and e = 0.5. 



Figure 8: (Color online) Bifurcation diagram of localized structures for model Eq. (1), in parameter plane (a, e), with cr = 0.8 
and n = 2. 


competition e, the uniformly vegetated state becomes unstable with respect to symmetry breaking instability. The 
threshold associated with this instability is represented by a solid line. This line is obtained by plotting formula Eq. 
(5). This spatial instability impedes the existence of fairy circles in the region C, and may allow for the formation of 
periodic structures. For a > 1/2, a single fairy grows until infinity and disappears, and a localized vegetation patch 
appears. This structure is stable in the region LVP, as shown in the bifurcation diagram of Fig. 8. The region B 
corresponds to a high degree of aridity. In this zone of parameters, a localized patch shrinks and disappears, and 
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Figure 9: (Color online) Two dimensional localized vegetation patterns obtained from numerical simulations of Eq. (1). (a) 
Multi-fairy circles for a = 0.44 (b) Multi-localized vegetation patches for a = 0.57. Other parameters are a = 0.1, e = 10 and 
n = 2.5. 


transition toward a bare state occnrs. 


V. CONCLUSIONS 

We have investigated the role of a strong nonlocal coupling in a bistable model namely the Nagumo model. This 
prototype model of population dynamics could be applied to vegetation dynamics. We have shown that far from any 
symmetry breaking or Turing type of instability, localized vegetation structures can be stabilized in large values of 
the aridity parameter. Their formation is attributed to the interaction between fronts mediated by a strong nonlocal 
coupling in the form of a Lorentzian. We have identified the following scenario: when increasing the level of the aridity, 
a large dip embedded in a uniformly vegetated state is formed. This structure has a single fringe peak that appeared 
in the spatial profile of the biomass. We have interpreted this behavior as a fairy circle. When increasing further 
the degree of aridity, localized vegetation patch can be formed in the system. This structure has a peak surrounded 
by the bare state. The localized structures reported in this work have a varying width as a function of aridity. In 
contrast, the width of localized vegetation structures found close to the symmetry breaking instability is determined 
by the most unstable Turing wavelength [26, 70]. We have established analytically a formula for the width of fairy 
circles and localized vegetation patches as a function the degree of aridity. The width of these localized structures is 
intrinsic to the dynamics of arid ecosystem and it is independent of external environmental effects, such as termites 
or ants. The results of direct numerical simulations of model Eq. (1) agreed with the analytical findings. 

In this paper we have focused our analysis on a single localized structure, several of them could be stable as shown 
in the Fig. 9. The formation of multi-dips or peaks localized structures, their interactions and their stability are 
under investigation. Understanding the formation of localized structures is central not only in arid ecosystems but 
also in spatially extended out of equilibrium systems. 
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